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4. COMPOSITE THREE-DIMENSIONAL PROFILES OF THE VELOCITY AND PRESSURE 

FIELD IN THE A.F.C. ROTOR 

The flow field in the AFC rotor can be divided into different 
regions (inviscid core region, end-wall region, tip- leakage, blade 
boundary layer) . Different types of probes were used to measure the 
various regions. The type of probe used at each region is dictated 
by the flow field characteristics of that region (three or two 
dimensional, high or low turbulent) . Four different sets of measure- 
ments were performed in the PSU/TURBO AFC rotor. The core inviscid 
flow was measured with a rotating five hole probe, the end-wall 
region with a rotating three sensor hot wire, the blade boundary 
layer with a rotating two sensor hot wire and the tip leakage flow 
with a stationary two sensor hot wire. 

We are now in the process of combining all these data together 
in order to acquire a complete picture of the flow field in the rotor. 

Some of the completed composite profiles are shown in Figs. 16 through 20. 

Figure 16 shows the relative streamwise velocity plotted versus 
the tangential and radial directions at the streamwise location S = 0,979. 
This figure combines the blade boundary layer data, the inviscid core 
dcta and the end-wall data. In Fig. 17 the tip-leakage data are also 
included. As expected the relative velocity goes to one near the 
casing. 

Figure 18 shows the radial velocity (Q^) plotted versus the axial and 
tangential distance at R = 0.973. A strong inward radial velocity 
can be clearly seen at the mid-passage. This is because of the tip 
leakage. In this figure only the end-wall data are shown. 
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PREFACE 


The progress of research on "Three Dimensional Flow Field Inside a 
Compressor Rotor Blade Passage, Including Blade Boundary Layers" (NASA 
Grant NSG 3266) for the six-month period ending December 31, 1982, is 
briefly reported here. The effort on turbulence modelling is incorporated 
into the paper "A Turbulence Model for Three Dimensional Turbulent 
Shear Flow Over Curved Rotating Bodies," AIAA Paper No. 83-0359. A 
copy of this paper was transmitted to NASA recently. For the sake of brevity, 
the material in this paper will not be repeated here. 


B. Lakshminarayana 
Principal Investigator 
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nomenclature for CHAPTER 1 

Jacobian of the transformation 


r,0 ,z 
U,V,W 


n z ’ n e ,n r 


W R r 


^z’^G’^r 


C » n ,R 


Pressure 

cylindrical coordinate system 

velocities in axial, tangential, and radial direction, 
respectively 

3n 3n 3n 
3z’ r90 ’ 3r 

3R 3R_ 3R 
3z ’ r30 3r 

3£ 3£_ H 

3z* r30 3r 

molecular viscosity 

transformed coordinates in the streamwise, normal and 
radial directions, respectively w 

density 

angular velocity 
turbulent kinetic energy 
turbulent dissipation rate 
Reynolds number 
friction velocity 

resultant velocity parallel to the wall at the first grid 
point away from the wall 

internal energy 


u « n + u + kj 

Z U 1. 

V = n U + n c v + n W 
z o r 


w = R z U + RgV + R r W 


kinematic viscosity 


Subscripts 


value of the first grid point away from the wall 


V 


NOMENCLATURE FOR CHAPTERS 3 and 4 


C chord length 

£,d length and diameter of the hot wire 

N n/C 


R 

S 


r/r 


tip 


s/C 


PS,SS 


pressure and suction surface, respectively 


s.n.r 


U 

U 

I 

w 


T 

se 


T 

s 


T 

r 



streamwise, normal, and radial directions (orthogonal to 
each other) shown in Fig. 1. s = 0 at the leading edge, 
n • 0 on the blade surface, r ° 0 at the axis of the machine. 

streamwise relative velocity normalized by U g 

local free stream (or edge) relative velocity 

4 

the difference between the local radial velocity and the 
free stream radial velocity normalized by U . 

the free stream value of the streamwise relative intensity 

normalized by U . 

e 

the free stream value of the radial intensity normalized 

bv U . 
e 

relative streamwise intensity normalized by T 

se 

radial intensity normalized by T re 

relative streamwise velocity in the s-n plane normalized 
by the blade tip speed 


Qr 


Y 


radial velocity normalized by the blade tip speed 

tangential distance normalized by the spacing and measured 
from the blade surface (Y = 0 on the suction side, Y = 1 
on the pressure side). 

static pressure normalized by -|pU ^ 

total pressure normalized by -^pU ^ 


blade tip speed 


\\ 




I. NUMERICAL SOLUTION OF THE FLOW FIELD INSID' THE PASSAGE OF A 

TURBOMACHINERY ROTOR ?P ^AGE 

The space marching code, developed by Govindan and Lakshminarayana [1], 
has been modified in order to be able to predict the flow field inside a 
rotor passage, including the blade and hub wall boundary layers. The 
basic changes incorporated are as follows: 

(i) Modifications of the equations so that the code can handle 
three-dimensional configurations with changes in the radial 
direction (for example changes in stagger angle, blade camtar 
and thickness. * 

(ii) Extensions and modifications in order to implement physically 
realistic turbulence model su as k-e model and algebraic 
Reynolds stress model 

1 . 1 Modification of the Equations 

The equations in the original code are based on the transformation 
of the compressible flow equation from a cylindrical coordinate system 
(r,0,z) to a body fitted coordinate system (R,£,n) through the trans- 
formation 

i = £(z,e) 

n = n(z,0) (1) 

R = R(r) 

for the code to handle fully three dimensional geometries the transforma- 
tion should be 

K = E(r,0,z) 

n = n(r ,0 ,z) (2) 

R = R(r ,0 ,z) 


1 


2 


The equations and the corresponding Jacobian matrices have been 
derived in the body fitted coordinates given by Eq . (2). The results 
are given in Appendix A. 


1.2 Turbulence Modeling 

The algebraic Reynolds stress model developed in Ref. 2 and the 
k-e model described in Ref. 3 have been coded and integrated with the 
main space marching code. 

The calculation of the turbulent kinetic energy, dissipation and 
Reynolds stresses are lagged one streamwise step. The five equations 
(continuity, three momentum, and energy) are first solved using the 
turbulence quantities derived at the previous streamwise step. When 


p, U, V, W, e^ are calculated, the k and e equations are solved in 
order to derive the values of k and £. The algebraic system of 
equations for the Reynolds stresses is then invert 1 to derive the 
Reynolds stresses. These stresses are used in the next streamwise 
station . 

in order to avoid using a very large number of grid points near 
the walls, the so-called wall-functions are used [4]. This approach 
assumes that, at the first grid point away from the wall with wall 
distance n^ just outside the viscous sublayer, the velocity components 
parallel to the wall follow the logarithmic law of the wall and the 
turbulence is in local equilibrium. With these assumptions, the 
resultant velocity parallel to the wall Q^, the kinetic energy k^ 
and the dissipation rate at point n^ are related by the following 


relationships 


ORIGINAL PAGE S3 
OP POOR QUALITY 




in E 


where K, E, are constants; K = 0.4, E * 9, C = 0.09 

1.2.1 Status of the code The subroutines implementing the 
k-e model and the Reynolds stress model have been completed and 
integrated with the main code The code is presently tested in 
simple turbulent flow configurations for debugging purposes. 
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2. A TURBULENCE MODEL FOR THREE DIMENSIONAL TURBULENT SHEAR FLOW 

OVER CURVED ROTATING BODIES 

It is known that the curvature and rotation affect the turbulence 
structure substantially and a knowledge of these effects are essential 
for the improved prediction of flow over rotating bodies. A turbulence 
model which includes the effects of curvature as well as rotation was 
developed during this reporting period. For the sake of brevity, this 
analysis is not presented in this report as it has been published 
as an AIAA Paper, "A Turbulence Model for Three Dimensional Shear Flow 
Over Curved Rotating Bodies," J. Galmes and B. Lakshminarayana , AIAA 
Paper No. 83-0559. Different hypotheses introduced to model the higher 

j 

order unknowns in the Reynolds stress equations are presented in this 
paper; a set of algebraic equations is derived. The transport 
equations of the turbulent kinetic energy and dissipation rate are 
discussed. A detailed analysis of the effect of the rotation on 
each component of the Reynolds stress tensor is presented for 
hypothetical cases such as the pure shear flow in a rotating frame. 
Calculations show that the effects of rotation on turbulent shear 
stresses are more pronounced in a centrifugal type of r.urbomachinery 


than an axial type. 
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3. BLADE BOUNDARY LAYER IN AN AXIAL 
FLOW COMPRESSOR ROTOR 

Th-?. three-dimensional turbulent boundary layer developing on 
the rotor of the compressor of the PSU/TURBO Lab was measured using 
a miniature "x" configuration hot wire probe. This investigation 
started earlier and its major part was completed in the period 
July - December 1982. The measurements were carried out at nine 
radial locations on both surfaces of the blade at various chordwisc 
locations (see Table 1) . 

t 

3 . 1 Measurement Technique and Corrections for Wall Vicinity Effect 
All the velocity and turbulence measurements were taken with a 
miniature crossflow "x" wire probe TSI (1247) shown as an insert in 
Fig. 2. The sensors were 3 mm diameter platinum- tungsten wires with 
£/d % 300. The sensors were located in the (sr) plane with their 
axis at 45° to the s axis (Fig. 1) and were traversed normal to the 
blade surface. Since the flow traverse was done close to the blade 
surface, the component of velocity i n the n direction (V) is assumed to 
be small. The present configuration together with the relevant hot 
wire equations [5] provide the value of the streamwise velocity U, 
radial velocity W and the respective intensities T^ and T^ in the 
s,n,r system shown in Fig. 1. The s coordinate is parallel to the 
blade surface lying on the cylindrical plane, n is the principal normal 
and r is the radial direction as shown in Fig. 1. The coordinate 


system is orthogonal. 


Since the measurements were taken very close to the blade surface 


and the heat transfer characteristics are affected by the wall vicinity, 
the hot wire probe was calibrated To derive corrections for the wall 
vicinity effect. These are incorporated int.o the hot wire equations of 
Ref. 5. The probe was calibrated in a jet with a wall (parallel to the 
jet) at the exit. The distance between the wall and the probe contin- 
uously varied from 0.05 mm to 5 mm, beyond which the change in 
calibration curve was found to be negligibly small. The calibration 
curve is shown in Fig. 2. The data includes the voltage at zero 
velocity (E o ) as well as the voltage (E) at various jet velocities (Q) . 
It is clear from Fig. 2 that the exponent in King's law is not 
affected by the wall. King's law is given by i 

E 2 - E 2 = BQ n (3) 

o 

where B is a calibration constant and n is the exponent. The values of 

E and B are affected by the wall vicinity. The correlation for the 
o 

wall vicinity based on this data is given by. 


2 2 2 D 2 

E Z = (F - E Z ) + E Z . 

C Odo ow D ow 


(4) 


where E is the value of E at D = D , and E is the voltage at zero 
Oco o 00 ow ° 

velocity near the wall. D is the distance between the sensor nearest 
to the wall and the wall. The value of D = 5 mm and D = 0.05 mm, 

1.25 mm, 2.5 mm, and 5 mm. The shift in the data shows that the 
coefficient B is affected and the following correlation for B is 
derived from this data. 


3 = C B^ exp(- ~ log c) , C = 1.01 


where B is the coefficient away from the wall. 
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It appears that without the corrections introduced through Eqs. 
(4) and (5) there is an e'ror of approximately 6 percent in the 
velocity close to the wall. These corrections (Eqs. (4) and (5)) have 
been incorporated in the hot wire equations. The value of E qw and 
E were measured at the beginning and the end of each measurement 

Oco 

run. 


3.2 Correction for Spatial Error 

The two sensors of the probe are not exactly at the same location 
but are spaced at 0.51 mm apart. In order to correct the error 
introduced by this distance, the velocity that each wire feels 
(V 1 ,V 9 ) is plotted separately versus the true distance of each wire 
(N ,N 0 ) from the -Jail (Fig. 3). Then the streamwise and radial 
velocity components at some distance re distance of the closest 

to the wall wire) are calculated usinf, the value of at and the 
interpolation (or extrapolated 1 ) values of at the location . 

Figure 4 shows both the corrected and uncorrected values of U and W 
for the measured boundary layers on the suction side of the blade 
at R = 0.75. All the subsequent figures show rhe corrected values 
of U and W. The turbulence intensities cannot be corrected for 
the spatial error. Judging from the mean velocites (Fig. 4) the 
correction is not large and it is confined only to a few points near 
the wall. 

3 . 3 Experimental Results 

Results at selected radial locations, R = 0.75, 0.918, 0.98, are 
given. The complete set of data will be included in a report under 
preparation . 


3 


Figures 4 through 6 show the boundary layer development on the 
suction side of the blade at the three radial locations. Figures 7 
through 9 show the corresponding results on the pressure side of the 
blade. It can be clearly seen that at the trailing edge region the 
boundary layer thickness is approximately two times larger on the 
suction side than on the pressure side.. Going from lower to higher 
radii the boundary layer thickness is increasing on the suction side 
and stays approximately constant on the pressure side. In the end- 
wall region (R = 0.918 to R ° 0.98), the effect of tip leakage which 
acts as a "suction" on the boundary layer on the pressure side and 
as a "blowing" on the suction side can be teen. 

The radial velocity is outwards at most locations. On the 4 
suction side it is decreasing towards the top while on the pressure 
side it is increasing. This again can be attributed to the tip 
leakage effects. 

Figures 10 through 15 show the turbulent quantities at the 

three radial locations on the suction and pressure sides of the 

blade. The turbulent intensities T , T are normalized respectively 

by the free stream turbulent intensities T ge , T . It is interesting 

to note t ... t T and T have almost the same values at all locations, 
r s 

For a stationary boundary layer T g > T^ . The reason that T g ^ T^ 
is probably due to the rotation effect, analyzed in Ref. 6. 


10 

Figures 19 and 20 show, respectively, the static pressure P g 
and tht total pressure plotted versus the radial and tangential 
directions. Pg is at the strearawise location S “ 0.25 and P^, is at 
S « 0.979. In the last two figures only the inviscid core data are 
shown. 








Table 1 

Radial and Streamwise Measurement Locations 
of the Blade Boundary Layer 


R 

Pressure Side 
Distance S 

Suction Side 1 

Distance S j 

0.583 

0.22 

0.44 

0.68 

0.92 

0.49 

0.66 

0.81 

0.99 

0.67 

0.23 

0.48 

0.70 

0.93 

0.55 

0.69 

0.84 

0.99 

0.75 

0.24 

0.52 

0.73 

0.94 

0.60 

0.72 

0.87 

0.99 

0.823 

0.245 

0.525 

0.75 

0.95 

0.63 

0.77 

0.9 

t) .99 

0.918 

0.25 

0.53 

0.76 

0.96 

0.65 

0.81 

0.94 

0.99 

0.945 

0.26 

0.54 

0.78 

0.97 

0.66 

0.83 

0.96 


0.959 

0.265 

0.545 

0.785 

0.975 

0.665 

0.835 

0.97 


0.973 

0.27 

0.55 

0.79 

0.98 

0.67 

0.84 

0.98 


0.98 

0.27 

0.55 

0.79 

0.98 

0.67 

0.84 

0.99 






12 


PUBLICATIONS & PRESENTATIONS 


The following papers were published during this period: 


(1) P. Lakshminarayana, T. R. Govindan, and C. Hah, "Experimental 
Study of the Boundary Layer on a Turbomachinery Rotor Blade." 
In Three Dimensional Turbulent Boundary Layers , edited by 
H. H. Fernholz and E. Krause, IUTAM Symposium Proceedings; 
Springer-Verlag, pp. 165-176. 


The following presentations were made during the reporting 
period: 


(1) "Computation Measurement of Rotor Wake Including Turbulence 
Modelling," William Maxwell Reed Mechanical Engineering Seminar , 
University of Kentucky, September 23, 1982. 

(2) "Turbulence Modelling and Turbomachinerv Flow Computation," Workshop 
on Trubomachinery Flow Computation, NASA Lewis Reserach Center, 
October 21, 1982. 
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The continuity, the three momentum and the energy equations are 
written in the form 


31 + 12 + 2£ + c = -1 

35 + 3n 3R + C Re 


|| + f + sl 

3R 9n J 


E = - 


E, F, G, C, P, Q, S are 5x1 column vectors and are given below: 


pv 


pvll + n P 
z 


G = - 


e, : 

F, 

G, C, 


Pu 


puU 

+ 

V 

puV 

+ 

V 

puW 

+ 

5 r p 

pufy 

6 i 


- 

pw 


pwU 

+ 

R P 
z 

pwV 

+ 

R 0 P 

pw W 

+ 

R P 
r 

pw (y 

'e . 

l 



F = 4 


C = 


pvV + n^P 


pvW + n^P 


pv 


(Y e i 


2 2 2 
U +V +W • 


W/R 

uw 


2p 


vw 


R 


■+ 2 P nw 


2 2 

(W -V ) 2 

p , S_w — - L - Z . _ p n z R - 2pfiV 

R 

2 2 2 

W, N . U +V +W ^2„ TT 
-(ype ) + p pft RW 


Pi - 0 

p = k [ (R 2 + 2 4 r 2 |2 + || + |S + (R n +R 4 R )B 

2 J [ r 0 3 z 3R 3 3R 3 dR rr 003zz3n 


3 V 


3W 


W 


+ (R^n - ~R n + (R n -^Rn)^-^R^ 
bz 3z03n rz 3zr3n 3 z R 



safes; 
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3 J 


R e R z 3U . /n 2 , 4 „2 S 3V , R r R 0 3W , „ 2 „ x 3U 

— t — Tn + (R + "T Rq+ R ) irr + — = — Trr + (nR - — R n ) -r— 

3 3R r 3 0 z 3R 3 3R 0 z 3 0 z 3n 
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5 J 


IE 

3n 


/ n n 

(U(n R + n„R 0 + - n R ) + V(R n n - -r R n ) + W(R n - - R n )} 
v rr 00 3zz z 0 30z z r 3 r z ; 


+ |^(u(R 2 + R 2 + 4 R 2 ) + 4 R fi R v + i R R ^ + 3 v fv(R n + 4 Ron- + R n ) 
3R V r 0 3z 30z 3 r z - — ' r r 3 t) 0 z z 

3n 


+ u(V* - f R z n e ) + "<V r - ! + lf( v(E r + ! R e + R P ‘ 


R R, 


+ U + W) + n R + HqRq + n R ) + U(n R - 4 a R ) 

3 3 1 3n v 3 r r 0 o zz zr 3 r z 


R„R 


0 r 


3 V r,,A 


R R 


R q R_V 


+ V(n 0 R r ‘ ! E 0 n r^ + UK R r + E e 

+ -z—fCR n + n R + n R ) 
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i o VWR fl „ 

+ ^ ("V Z R + -^-2- - 4 UWR ) 
R r 3 3 z y 


^ e i 2 2 
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Qi - o 
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The Jacobian Matrix of the Vector E 
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